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Abstract 

If the background universe is observed through a significant amount of diffusely dis¬ 
tributed foreground dust, then studies at optical wavelengths may be severely biased. Pre¬ 
vious studies investigating the effects of foreground dust on background sources assumed 
dust to be “compactly” distributed, ie. on scales comparable to the visible extent of nor¬ 
mal galaxies. We show however that diffuse dust is more effective at obscuring background 
sources. Galaxy clusters are a likely location for “large-scale” diffusely distributed dust, and 
its effect on the counts of background sources is explored. We also explore the implications 
of a hypothesised diffuse intergalactic dust component uniformly distributed to high redshift 
with comoving mass density equal to that associated with local galaxies. In this case, we 
predict a deficit in background sources about three times greater than that found in previous 
studies. 
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1 Introduction 

There are a number of studies claiming that dust in foreground galaxies has a substantial effect on 
the colours and counts of optically-selected quasars (Ostriker & Heisler 1984; Heisler & Ostriker 
1988; Fall & Pei 1992 and Wright 1990). It is estimated that at least 50% of bright quasars 
at redshifts z>3 may be obscured by foreground galactic dust and hence missing from optical 
samples. These studies assumed that dust was conhned only within the visible extent of normal 
massive galaxies. However, distant populations such as faint field galaxies and quasars may also 
be observed through foreground diffuse dust distributions. Such distributions may be associated 
with galaxy clusters and extended galactic haloes. 

A truly diffuse, intergalactic dust component is ruled out based on the counts of quasars and 
reddening as a function of redshift (eg. Rudnicki 1986; Ostriker & Heisler 1984). Such observations 
indicate that if a signihcant amount of dust exists, it must be patchy and diffuse with relatively 
low optical depth so that quasars will appear reddened without being removed from flux-limited 
samples. 

Galaxy clusters provide a likely location for ‘large-scale’ diffusely distributed dust. Indirect 
evidence is provided by several studies reporting large dehcits of distant quasars or clusters of 
galaxies behind nearby clusters (Boyle et al. 1988; Romani & Maoz 1992 and references therein). 
These studies proposed that extinction by intracluster dust was the major cause. Additional 
evidence for diffuse dust distributions is provided by observations of massive local galaxies where 
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in a few cases, dust haloes extending to scales >50kpc have been conhrmed (Zaritsky 1994 and 
Peletier et ah 1995). 

Does uniformly distributed dust really exist in the intergalactic medium (IGM)? Galactic winds 
associated with prodigious star formation at early epochs may have provided a likely source of 
metal enrichment and hence dust for the IGM (eg. Nath & Trentham 1997). Observations of metal 
lines in Ly-a absorption systems of low column density (iVni^lO^^cm”^) indeed suggest that the 
IGM was enriched to about Z ~ O.OIZq by redshift z ~ 3 (Womble et ah 1996; Songaila & Gowie 
1996). A source of diffuse dust may also have been provided by an early generation of pre-galactic 
stars (ie. population III stars) associated with the formation of galactic haloes (McDowell 1986). 
Studies have shown that possible reddening from uniformly distributed IGM dust is limited by 
observations of radio-selected quasars. Since radio-selected quasars should have no bias against 
reddening by dust, such a component must be of sufficiently low optical depth to avoid producing 
a large fraction of ‘reddened’ sources at high redshift (see Webster et ah 1995; Masci 1997). 

In this paper, we show that a given quantity of dust has a much greater effect on the back¬ 
ground universe when diffusely distributed. We shall investigate the effects of diffuse dust from, 
hrst, possible distributions on galaxy cluster scales and second, from a hypothesised uniformly 
distributed component in the IGM. 

This paper is organised as follows: in the next section, we explore the dependence of background 
source counts observed through a given mass of dust on its spatial extent. In Section |^, we 
investigate the spatial distribution of dust optical depth through galaxy clusters and its effect 
on the counts and colours of background sources. Section || explores the consequences if all dust 
in the local universe were assumed uniformly distributed in the IGM. Further implications are 
discussed in Section ^ and all results are summarised in Section |^. All calculations use a Friedmann 
cosmology with qq = 0.5 and Hubble parameter h .50 = 1 where Hq = hOhsokms"^ Mpc“^. 

2 Compact versus Diffuse Dust Distributions 

In this section, we explore the dependence of obscuration of background sources on the spatial 
distribution of a given mass of dust. For simplicity, we assume the dust is associated with a 
cylindrical face-on disk with uniform dust mass density. We quantify the amount of obscuration 
by investigating the number of background sources behind our absorber that are missed from an 
optical flux-limited sample. 

The fraction of sources missing to some luminosity L relative to the case where there is no dust 
extinction is simply l — Nobs{> L)/Ntrue{> L), where Nobsi> represents the observed number of 
sources in the presence of dust. For a uniform dust optical depth r, Nobs{> L) = Ntrue{> e^L). For 
simplicity, we assume that background sources are described by a cumulative luminosity function 
that follows a power-law: 4) oc where (5 is the slope. This form for the luminosity function 
is often observed for ‘luminous’ (L > L*) galaxies and quasars, which dominate high redshift 
( 2 ;> 0 . 5 ) populations in flux-limited samples. With this assumption, the fraction of background 
sources missing over a given area when viewed through our dusty absorber with uniform optical 
depth r is given by 

fmissir) = 1 - (1) 

If the ‘true’ number of background sources per unit solid angle is ntme, then the total number of 
background sources lost from a flux limited sample within the projected radius R of our absorber 
can be written: ^ 

Nlostir) = fibrue'^il - e-^n, (2) 
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where D is the distance of the absorber from us. 

To investigate the dependence of background source counts on the spatial dust distribution, 
we need to hrst determine the dependence of r in equation on the spatial extent R for a fixed 
mass of dust Md- This can be determined from the individual properties of grains as follows. The 
extinction optical depth at a wavelength A through a slab of dust composed of grains with uniform 
radius a is dehned as 

Tx = Qext{Ka)7ia‘^ndld, (3) 

where Qext is the extinction efficiency which depends on the grain size and dielectric properties, 
Ud is the number density of grains and Id is the length of the dust column along the line-of-sight. 
Assuming our cylindrical absorber (whose axis lies along the line-of-sight) has a uniform dust mass 
density: pdust = Md/nR'^ld, where R is its cross-sectional radius, we can write, Ud = Pdust/^'^o.^Pd, 
where pd is the mass density of an individual grain. We use the extinction efficiency Qext ffi the 
Id-band as parameterised by Goudfrooij et ah (1994) for a graphite and silicate mixture (of equal 
abundances) with mean grain size a ~ O.l/rm characteristic of the galactic ISM. The value used 
is Qext{V, 0.1pm) = 1.4. We use a galactic extinction curve to convert to a i?-band extinction 
measure, where typically tb ~ 1.3ry (eg. Pei 1992). Combining these quantities, we hud that the 
S-band optical depth, tb, through our model absorber can be written in terms of its dust mass 
and cross-sectional radius as follows: 

Md \ f R r' f a ( Pd 
IO^MqJ y20kpcy yO.l/rmy y2gmcm“^ 

where we have scaled to a dust mass and radius typical of local massive spirals and ellipticals (eg. 
Zaritsky 1994). This measure is consistent with mean optical depths derived by other means (eg. 
Giovanelli et al. 1994 and references therein). 

From equation (Q), we see that the dust optical depth through our model absorber for a fixed 
dust mass varies in terms of its cross-sectional radius i? as r oc 1/R^. For the nominal dust 
parameters in equation (H), the number of sources missed behind our model absorber (equation 
^ can be written 

Niost{< R) = Ntrue{< 20kpc) RIq{1 - (5) 

where i? 2 o = i?/20kpc and 

/ IN 7r(20kpc)^ , , 

Ntrue{< 20kpc) = fltrue -- (6) 

is the ‘true’ number of background sources falling within the projected scale radius R = 20kpc. 
From the functional forms of equations (|]) and (j^), there are two limiting cases: 

1. For optical depths r 3> 1//?, the factor 1 — in equation (||) is of order unity. This 

corresponds to values of R such that 0 < i? 2 o for the nominal parameters in equation 

(I). For values of R in this range, we have Wost oc R^ and the obscuration of background 
sources will depend most strongly on R. 

2. For r -C 1//3 or equivalently i? 2 o S> Niost will approach a constant limiting value, 

independent of the dust extent R. From equation (^), this limiting value can be shown to 
be Niosti< R) = Ntruei< 20kpc)/3. 

As a simple illustration, we show in Figure |l| the dependence of the number of background 
sources missing behind our model dust absorber on R, for a hxed dust mass of IO^Mq as dehned 
by equation (||). We assumed a cumulative luminosity function slope of = 2.5, typical of that 
for luminous galaxies and quasars. From the above discussion, we see that when ie. 
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Figure 1: The number of sources missing from an optical flux-limited sample behind a face-on 
dusty disk with fixed mass M^ust = IO^Mq as a function of its radial extent R. This number is 
normalised against the ‘true’ number of background sources (in the absence of the absorber) that 
fall within the solid angle 7r(20kpc)^/Zi)^ subtended by our nominal radius of 20kpc at distance D. 
See equation (|[) 


i?>30kpc (or when = 0.4), the obscuration will start to approach its maximum value and 

remain approximately constant as i? —> cxo. 

We conclude that when dust becomes diffuse and extended on a scale such that the mean 
optical depth r through the distribution satishes r < 1//?, where f3 is the cumulative luminosity 
function slope of background sources, obscuration will start to be important and is maximised 
for r -C 1//9. The characteristic spatial scale at which this occurs will depend on the dust mass 
through equation (^. For the typical grain values in equation (|^), this characteristic radius is 
given by 


R ~ 



Md 

IQSMo 
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The simple model in Figure shows that the obscuration of background sources due to a normal 
foreground galaxy will be most effective if dust is distributed over a region a few times its optical 
radius. This prediction may be difficult to conhrm observationally due to possible contamination 
from light in the galactic absorber. In the following sections, we explore two examples of possible 
‘large-scale’ diffuse dust distributions that can be explored observationally. 


3 Diffuse Dust in Galaxy Clusters 

A number of studies have attributed the existence of large dehcits of background sources behind 
nearby galaxy clusters as due to extinction by dust. Bogart & Wagner (1973) found that distant 
rich Abell clusters were anticorrelated on the sky with nearby ones. They argued for a mean 
extinction of Ay k, 0.4 mag extending to ~ 2.5 times the optical radii of the nearby clusters. 
Boyle et al. (1988) however claimed a ~ 30% dehcit of background quasars within 4' of clusters 
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consisting of tens of galaxies. These anthors attribnte this to an extinction Ay ^ 0.15 mag, and 
dednce a dust mass of ~ 10 ^*^Mq within 0.5Mpc of the clusters. Romani & Maoz (1992) found that 
optically-selected quasars from the Veron-Cetty & Veron (1989) catalogue avoid rich foreground 
Abell clusters. They also found dehcits of ~ 30% out to radii ~ 5' from the clusters, and postulate 
a mean extinction. Ay ^ 0.4 mag. 

The numbers of background sources behind clusters is also expected to be modihed by gravi¬ 
tational lensing (GL) by the cluster potential. Depending on the intrinsic luminosity function of 
the background population, and the limiting magnitude to which the sources are detected, GL can 
cause either an enhancement or a dehcit in the number of background sources. The GL effect has 
been used to explain various reports of overdensities of both optically and radio-selected quasars 
behind foreground clusters (Bartelmann & Schneider 1993; Bartelmann et al. 1994; Rodrigues- 
Williams & Hogan 1994; Seitz & Schneider 1995). The reported overdensities for optically-selected 
QSOs are contrary to the studies above where anticorrelations with foreground clusters are found. 
These overdensities however are claimed on angular scales ~ 10' — 30' from the cluster centers, 
considerably larger than scales on which most of the under densities have been claimed, which are 
of order a few arcminutes. One interpretation is that dust obscuration bias may be greater towards 
cluster centers due to the presence of greater quantities of dust. On the other hand, the reported 
anticorrelations on small angular scales can perhaps be explained by the optically crowded helds, 
where QSO identihcation may be difficult. At present, the effects of clusters on background source 
counts still remains controversial. 

More direct evidence for the existence of intracluster dust was provided by Hu et al. (1985) 
and Hu (1992) who compared the Ly-o; flux from emission line systems in “cooling flow” clusters 
with Balmer line fluxes at optical wavelengths. Extinctions of Ay ~ 0.2 — 1 mag towards the 
cluster centers were found, in good agreement with estimates from the quasar dehcits above. 
Intracluster dust has been predicted to give rise to detectable diffuse IR emission (eg. Dwek et 
al. 1990). The most extensive search was conducted by Wise et al. (1993), who claimed to have 
detected excess diffuse 60-100/im emission at the 2a level from a number of rich Abell clusters. 
They derived dust temperatures in the range 24-34 K and dust masses ~ IO^^Mq within radii 
of ~ IMpc. Recently, Allen (1995) detected strong X-ray absorption and optical reddening in 
ellipticals situated at the centers of rich cooling how clusters, providing strong evidence for dust. 
These studies indicate that intracluster dust is certainly present, however, the magnitude of its 
ehect on producing background source dehcits remains a controversial issue. 

In this section, we give some predictions that may be used to further constrain cluster dust 
properties, or help determine the dominant mechanism (ie. GL or extinction) by which clusters 
ahect background observations. 

3.1 Spatial Distribution of Cluster Dust 

X-ray spectral measurements show the presence of hot, metal enriched gas in rich galaxy clusters 
with ~ 0.5 — 1 solar metallicity. This gas is believed to be both of galactic and primordial origin (ie. 
pre-existing IGM gas), with the bulk of metals being ejected from cluster galaxies (see Sarazin 
1986 for a review). Ejection from galaxies may occur abruptly through collisions between the 
cluster galaxies, ‘sudden’ ram-pressure ablation, or through continuous ram-pressure stripping by 
intracluster gas (eg. Takeda et al. 1984). The lack of signihcant amounts of dust (relative to 
what should have been produced by stellar evolution) and interstellar gas in cluster ellipticals 
provides evidence for a mass loss process. On the other hand, in ellipticals that avoid dense 
cluster environments, signihcant quantities of neutral hydrogen, molecular gas and dust have been 
detected (eg. Lees et al. 1991). 
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If the dust-to-gas ratio of intracluster gas in rich clusters were similar to that of the Milky 
Way, then a radial gas column density of typically >10^^cm“^ with metallicity Z = 0.5Zq would 
produce an extinction v4y>4 mag. This however is much greater than the value observed. The 
likely reason for the dehciency of dust in the intracluster medium is its destruction by thermal 
sputtering in the hot gas, a process which operates on timescales tsput 4 < 10 ®n_ 3 ao.i yr, where 
n _3 = hh is the gas density and ao.i = a/O.l/rm the grain radius (Draine & Salpeter 
1979). Dust injection timescales from galaxies is typically of order a Hubble time (eg. Takeda 
et ah 1984) and hence, grains are effectively destroyed, with only the most recently injected still 
surviving and providing possibly some measurable extinction. 

The spatial distribution in dust mass density remains a major uncertainty. A number of authors 
have shown that under a steady state of continuous injection from cluster galaxies, destruction by 
thermal sputtering at a constant rate, and assuming instantaneous mixing with the hot gas, the 
resulting mass density in dust will be of order 

~ (oi^) (m) 

(eg. Dwek et al. 1990) where Zd is the injected dust-to-gas mass ratio, assumed to be equal to the 
mean value of the galactic ISM, Zd — 0.01 (Pei 1992). According to this simple model, the dust 
mass density is independent of gas density and position in the cluster. If we relax the assumption 
of instantaneous mixing of dust with the hot gas however, so that the spatial distribution of gas 
is different from that of the injected dust, the radial distribution of dust can signihcantly differ 
from uniformity throughout a cluster. Such a non-uniform spatial dust distribution may be found 
in clusters exhibiting cooling flows. If, as suggested by Fabian et al. (1991), most of the cooled 
gas remains cold and becomes molecular in cluster cores, then a relatively large amount of dust 
may also form, resulting in a dust distribution which peaks within the central regions. 

We explore the radial dependence of extinction optical depth through a cluster, and the ex¬ 
pected dehcit in background sources by assuming that dust is diffusely distributed and follows a 
spatial density distribution: 


Pdust{,R^ P(0) 



(9) 


where Rc is a characteristic radius which we £x and n is our free parameter. Equation (^) with 
u = 3/2 is the usual King prohle which with — 0.25Mpc, represents a good approximation to 
the galaxy distribution in clusters. Thus for simplicity, we keep Rc fixed at Rc = 0.25Mpc and 
vary n. To bracket the range of possibilities in the distribution of intracluster dust, we consider 
the range 0 < n < 3/2. n = 0 corresponds to the simple case where Pdust{R) = p(0) =constant, 
which may describe a situation where injection of dust is balanced by its destruction by hot gas 
as discussed above. The value n = 3/2 assumes that dust follows the galaxy distribution. This 
prohle may arise if grain destruction by a similar distribution of hot gas were entirely absent. 


3.2 Spatial Distribution of Dust Optical Depth and Background Source 
Deficits 

To model the spatial distribution of optical depth, we assume that intracluster dust is distributed 
within a sphere of radius Rmax- The central dust mass density p(0) in equation (^ is hxed by 
assuming values for the total dust mass Mdust and Rmax such that 

r Rmax _ 

Mdust = / 4:TtR^p{R) dR. (10) 

JO 
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We assume a cluster dust radius of Rmax = IMpc, which represents a radius containing >90% of the 
virial mass of a typical dense cluster characterised by galactic velocity dispersion ay ~ 800kms~^ 
(Sarazin 1986 and references therein). We assume that the total dust mass within Rmax = IMpc 
is Mdust = This value is consistent with that derived from extinction measures by Hu, 

Cowie & Wang (1985), IR emission detections by Wise et ah (1993) and theoretical estimates of 
the mean intracluster dust density as given by equation (H). 

Using equation (^, the R-band optical depth through our spherical intracluster dust distribu¬ 
tion at some projected distance r from its center can be written: 




3Q 


ext 


^apd 



Pdust {[r^ + dR', 
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where pdust{R) is our assumed radial density distribution (equation |) and pd the mass density 
of an individual dust grain. For a uniform dust density {pdust{R) = Pdust{0) =constant), and our 
assumed values of Rmax and Mdust given above, the radial dependence in dust optical depth can 
be written: 


TB{r) 


Tb{ 0) 



Rn 


1/2 


( 12 ) 


where 'rB(O) is the optical depth through the center of our cluster, which with grain properties 
characteristic of the galactic ISM, will scale as 


^b(O) 




0.06 


/ Mdust \ 


f Rn.ax \ V g \ V Pd \ ' 

\ Mpc J \ 0 . 1 /imy y 2 gmcm“^y 


(13) 


This value is about a factor three times lower than estimates of the mean extinction derived from 
the dehcit of QSOs behind foreground clusters (eg. Boyle et al. 1988), and that implied by the 
Balmer decrements of Hu et al. (1985). For a hxed dust mass of IO^^Mq however, we can achieve 
larger values for the central optical depth by steepening the radial dust-density distribution prohle, 
determined by the slope n in equation (^. 

Figure ^ shows the optical depth as a function of projected cluster radius for the cases n = 0, 
0.5, 1, and 1.5. The case n = 0.5 approximately corresponds to the model of Dwek et al. (1990), 
which included effects of mild sputtering by hot gas in order to £t for the observed IR emission 
from the Coma cluster. As shown, the case n = 0 {pdust{R) =constant) predicts that the dust 
optical depth should be almost independent of projected radius r. Within all projected radii, the 
optical depths predicted by our diffuse dust model lie in the range 0 < rs < 0.3. Turning back to 
the discussion of Section ^ where we show that background obscuration by diffuse dust reaches its 
maximum for < 1//?, these optical depths satisfy this condition for /?<2.5, typical of luminous 
background galaxies and quasars. 

We now explore the effects of these models on background source counts as a function of 
projected cluster radius. We hrst give an estimate of the projected radius at which the numbers 
of background sources lost from a flux-limited sample is expected to be a maximum. This is 
determined by investigating the dependence in the differential number of sources missing, dNiost, 
within an interval {r,r + dr) as a function of projected radius r. From equation (j^), this differential 
number will scale as 


dNiost(r) (X r dr 


1 


_ g-PRr) 


(14) 


where r(r) is given by equation (^). Figure plots dNiost/dr as a function of r for our various 
models, where we have assumed (3 = 2.5. Thus from observations, an identihcation of the pro¬ 
jected radius at which the background source dehcit peaks can be used to constrain the spatial 
distribution of intracluster dust. 
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Figure 2: (a) Optical depth as a function of projected cluster radius for various dust density 
distributions as parameterised by equation (^. (b) Differential number of background sources 
lost from a flux-limited sample (arbitrary scale), (c) Total fraction of background sources missing 
to some r and (d) Cluster-QSO two-point angular correlation function, filled circles: Boyle et 
al. (1988), squares: Romain & Maoz (1992), open circles: Rodrigues-Williams & Hogan (1994), 
triangles: Rodrigues-Williams & Hawkins (1995). 



































































The cumulative fraction of background sources missing within a projected cluster radius is 
given by 

Niost{< r) _ 1 

~ „2 




2Tir' [l - dr' 


irr^ 


(15) 


d^truei^^ 

This fraction is shown in Figure |^c. As expected, the n = 3/2 model which contains the largest 
amount of dust within the inner few hundred kiloparsecs predicts the strongest trend with r, while 
the opposite is predicted if the dust density were completely uniform. These predictions can be 
compared with a number of existing studies of the observed two-point angular correlation function 
between clusters and optically-selected QSOs. This function is usually dehned as 


w. 


cq 


{0) = 


{NoU< 0)) 

{Nrani< 9)) 


- 1 , 


(16) 


where {Nobs{< d)) is the average number of observed cluster-QSO pairs within an angular radius 9 
and {Nran{.< 9)) is that expected in a random distribution. For our purposes, {Nran{.< 9)) can be 
replaced by the “true” number of cluster-QSO pairs expected in the absence of dust, and hence, 
we can re-write equation ([I6| ) as 

/n\ _ (A^/ost(^ ^)) _ » / /I 


We compare our models with a number of studies of Wcq{9) for optically selected QSOs in 
Figure These studies considerably differ from each other in the selection of the QSO and 
cluster samples, and as seen, both anticorrelations and correlations on different angular scales are 
found. The former have been interpreted in terms of extinction by intracluster dust, while the 
latter with the GL phenomenon. In most cases, the reported overdensities are too large to be 
consistent with GL models given our current knowledge of cluster masses and QSO distributions. 

It is interesting to note that the studies which have reached the smallest angular scales (^5') are 
also those in which anticorrelations between QSOs and foreground clusters have been reported. 
This can be understood in terms of a larger dust concentration and hence extinction towards 
cluster centers. These studies however may not be free of selection effects, such as in the detection 
of QSOs from the visual inspection of objective prism plates. From a cross-correlation analysis 
of galactic stars with their cluster sample however, Boyle et ah (1988) found that such selection 
effects are minimal. 

The maximum dust radial extent assumed in our models, Rmax = IMpc, corresponds to 
angular scales ~ 5' at the mean redshift of the clusters {{zc) ~ 0.15) used in these studies. Thus, 
as shown in Figure |^d, our model predictions only extend to ~ 5'. As shown in this hgure, the 
n = 3/2 model which corresponds to the case where the dust density is assumed to follow the 
galaxy distribution, provides the best fit to the Boyle et al. (1988) data. We must note that this is 
the only existing study performed to angular scales ~ 1' with which we can compare our models. 
Further studies to such scales are necessary to conhrm the Boyle et al. result, and/or provide a 
handle on any selection effects. 


3.3 Summary 

To summarise, we have shown that for a plausible value of the dust mass in a typical rich galaxy 
cluster, obscuration of background sources will be most effective if dust is diffusely distributed on 
scales ~ IMpc. This conclusion is based on our predicted optical depths {tb < 0.3) satisfying 
our condition for ‘maximum’ obscuration: tb < 1//3 (see Section where typically /3<2.5 for 
luminous background galaxies and QSOs. 
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We have explored the spatial distribution in dust optical depth and background source dehcits 
expected through a typical rich cluster by assuming different radial dust density profiles. These 
predictions can be used to constrain cluster dust properties. A dust density distribution with 
u = 3/2 (equation |^) appears to best satisfy the ‘small scale’ cluster-QSO angular correlation 
study of Boyle et ah (1988). 

4 Diffuse Intergalactic Dust? 

There have been a number of studies claiming that the bulk of metals in the local universe 
had already formed by z ~ 1 (eg. Lilly & Cowie 1987; White & Frenk 1991; Pei & Fall 1995). 
Similarly, models of dust evolution in the galaxy show that the bulk of its dust content was formed 
in the first few billion years (Wang 1991). These studies suggest that the global star formation 
rate peaked at epochs z>2 when the bulk of galaxies were believed to have formed. Supernova- 
driven winds at early epochs may thus have provided an effective mechanism by which chemically 
enriched material and dust were dispersed into the IGM. As modelled by Babul & Rees (1992), 
such a mechanism is postulated to be crucial in the evolution of the ‘faint blue’ galaxy population 
observed to magnitudes B ~ 28. Nath & Trentham (1997) also show that this mechanism could 
explain the recent detection of metallicities Z ~ O.OIZq in low density Ly-a absorption systems 
at 2 ; ~ 3. Another source of diffuse IGM dust may have been provided by an epoch of population 
III star formation associated with the formation of galactic haloes (eg. McDowell 1986). 

What are the effects expected on background sources if all dust formed to the present day 
was completely uniform and diffuse throughout the IGM? In this section, we show that such 
a component will have a low optical depth and have an insignificant effect on the colours of 
background sources, but will be high enough to significantly bias their number counts in the 
optical. 

4.1 Comoving Dust Mass Density 

To explore the effects of a diffuse intergalactic dust component, we need to assume a value for 
the mean mass density in dust in the local universe. This density must not exceed the total mass 
density in heavy metals at the present epoch. An upper bound for the local mass density in metals 
(hence dust) can be derived from the assumption that the mean metallicity of the local universe 
is typically: Z ~ ^metais/^gas ~ 0.01 (ie. the ratio of elements heavier than helium to total 
gas mass), as found from galactic chemical evolution models (eg. Tinsley 1976) and abundance 
observations (Grevesse & Anders 1991). Gombining this with the upper bound in the baryon 
density predicted from big-bang nucleosynthesis (Olive et ah 1990) where figas^^baryon < 0.06/i^o^, 
it is apparent that 

^metalsiz = 0) < 6 X (18) 

Let us now compute the total mass density in dust used in previous studies that modelled the 
effects of dust in individual galaxies on background quasars. Both Heisler & Ostriker (1988) and 
Fall & Pei (1993) modelled these effects by assuming that dust in each galaxy was distributed 
as an exponential disk with scale radius tq ~ 30kpc and central face-on optical depth, tb = 0.5. 
The comoving mean mass density in dust (relative to the critical density) in these studies, given 
a comoving galaxy number density Uq = 0 . 002 h 5 gMpc~^, can be shown to be 

^dusto ^ 7.3 X IQ-^ho ( - — - ) (19) 

V0-002Mpc“V V^Okpc/ V0.5y ^ ^ 
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(see Masci 1997). This is consistent with the constraint in equation (|T8[) . Thus, as a working 
measure, we assume the comoving mass density defined by equation (|^ in the calculation that 
follows. 


4.2 Obscuration by Diffuse Intergalactic Dust 


If the dust mass density given by equation 0) is assumed uniformly distributed and constant on 
comoving scales to some redshift, the i?-band optical depth through a dust sheet of width dZ at 
redshift 2 ; in an observer’s frame can be written (see equation ^ 

, / \ -1 / \ -1 


(Itb — 7.2 X 10 h 


50 


n 


dustO 


pd 


7.3 X 10-V lOT/rm^ 


2gm cm 


-3 


dZ ' 
Mpc 


(l + z), (20) 


where for simplicity, we have assumed dust properties characteristic of the galactic ISM. The factor 
(1 + z) is due to our assumption of a 1/A dependence for the dust extinction law. This arises from 
the fact that light received in the 5-band corresponds to light of wavelength Xb/{)- + z) at redshift 
2 ;, which consequently suffers greater extinction. Although a 1/A law is not fully representative 
of that observed in the galactic ISM which includes the strong 2200A feature, this is a good on 
average ralation for the dust laws in many external galaxies (eg. Jansen et al. 1994). Such an 
assumption greatly simplifies the redshift dependence of extinction in an observer’s frame. With 
dl/dz = 6000Zi5o^(l -I- ^)“®/^Mpc (for a go = 0.5 and A = 0 cosmology), the total mean optical 
depth to some redshift in an observer’s 5-band will scale as 


tb{z) ~ 0.1 


a 


-1 


dust{) 


pd 


7.3 X 10-V lOT/rm^ 


2gm cm 


-3 


(1 + ^) 


- 1 / 2 ' 


( 21 ) 


This represents the total optical depth if all dust in the intervening galaxy model of Heisler & 
Ostriker (1988) were assumed uniformly distributed throughout the universe. 

Assuming dust is uniformly distributed to z = 2, the observed 5-band optical depth from 
equation (pT]) will be of order 

rs(U)(z = 2) ~ 0.04. (22) 

Using a galactic extinction law (Pei 1992), this corresponds to an extinction in 5 — 5 colour 
of Eb-r ~ 0.02mag. Thus, if background faint field galaxies and QSOs are observed through a 
uniform intergalactic dust distribution, their observed colours are not expected to be significantly 
affected. We now show however that the numbers of sources missing at such redshifts could be 
significantly greater than that claimed by previous studies which assume all dust to be associated 
with massive galaxies alone. 

If dust to some distance D covers an area of sky A and hence, has covering factor Cd = A/47r5^, 
the number of background sources lost from a flux-limited sample can be estimated from equation 
(1^). In general, the number of background sources at some redshift lost from an area of sky with 
dust covering factor Cd will scale as: 


E^lostiC) ^') EdfmissiC) z)) 


(23) 


where fmiss = 1 — is the fraction of sources missing per unit area. For a completely 

uniform dust distribution, C'rf(U) = 1, and to redshift 2 ; = 2, fmiss — 10% for [3 = 2.5. 

If dust were confined to individual galaxies along the line-of-sight however, their covering 
factor, assuming they follow a Poisson distribution is typically Cd — A^zCxp {—Nz), where Nz is 
the mean number of absorber intersections to redshift z: 


Nz ~ O.OlZigo 


no 


0.002Mpc 


-3 



(IT^)'-"^-! 


(24) 
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(see Heisler & Ostriker 1988). We have scaled to the nominal parameters assumed in the interven¬ 
ing galaxy model of Heisler & Ostriker (1988) (hereafter HO). In this model, we hnd a covering 
factor of only Od(HO) ~ 0.04 to z = 2. We can estimate the mean effective optical depth observed 
in the H-band through an individual absorber to ^ ~ 2 in the HO model by using the formalism 
of Section |^. For a fixed mass of dust, equation (^ implies that the product of the area (or cov¬ 
ering factor) and optical depth of a dust distribution: r x Cd, is a constant, depending on grain 
properties and dust mass alone. Using this relation, the observed effective absorber optical depth 
to z ~ 2 in the HO model, rB(HO), can be estimated by scaling from our values of and 

Orf(U) above for uniformly distributed dust: 

tb(HO) ~ rB(U)C-,(U)/C,(HO) = = i. (25) 

Using this value, the fraction of background sources missed by obscuration from an individual 
absorber is ‘effectively’ fmiss = 1 — exp(—2.5 x 1) ~ 91%. Combining these results, we hnd using 
equation (|23|) that the number of sources missing at z ~ 2 due to a uniform foreground dust 
distribution to be greater by a factor of q () 4 xV 9 i ~ 3 than that predicted by Heisler & Ostriker 
(1988). 

We must note that this estimate makes no allowance for possible evolution in dust content. 
Effects of foreground diffuse dust on source counts at 2 ; > 2 may be signihcantly reduced if 
appreciable evolution has occured. Effects of models where the dust content evolves have been 
explored by Masci & Webster (1998) 

4.3 Summary 

We conclude that the existence of a signihcant amount of diffusely distributed dust (eg. with mass 
density on comoving scales of order that observed in local galaxies) can enhance the number of 
background sources missing in optical samples. Due to its relatively large covering factor, diffuse 
dust predicts a reduction in optical counts at 2 : > 2 about three times greater than that claimed 
by previous studies. 

The colours of background sources are not expected to be signihcantly ahected. This implies 
that the use of background populations to measure a dihuse IGM dust component will be extremely 
difficult. 


5 Discussion 

In this section, we discuss some further uncertainties and implications regarding the existence of 
dihuse dust in the universe. 

First, the ehects of intracluster dust on background sources critically depends on the amount 
of dust present, and its spatial distribution. Regardless of the mechanism by which grains are 
injected into the intracluster medium from galaxies, it is possible that a signihcant fraction are 
destroyed in the injection process. Signihcant amounts of hot gas are also believed to exist in 
the ISM of cluster ellipticals (eg. Forman et al. 1985). This gas is expected to destroy grains 
on timescales <10®yr (Draine & Salpeter 1979), much shorter than injection timescales. Such 
destruction mechanisms can thus prevent the formation of signihcant quantities of dust. 

It is possible that the spatial distribution of intracluster dust is not ‘dihuse’ and uniformly 
distributed, but inhomogeneous. For example, Fabian et al. (1991) propose that if most of the 
cooled gas resulting from cluster cooling hows remains cold and becomes molecular, then this may 
provide suitable conditions for large amounts of dust to form. A clumpy dust distribution that 
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follows cooling flow filaments may result, reducing the effective dust covering factor and hence 
background source dehcit. These issues need to be addressed before attributing such dehcits to 
extinction by dust. 

The existence of a smooth IGM dust component also remains a major uncertainty. Due to their 
deep gravitational potential, Margolis & Schramm (1977) showed that it is unlikely for supernovae- 
driven winds to expel significant quantities of dust from a massive galaxy to large scales. For low 
mass galaxies however (eg. dwarfs). Babul & Rees (1992) show that this mechanism can be 
effective. Such systems are believed to comprise a majority of the ‘faint-blue population’ which 
show an excess ~ 20 — 30 times that predicted from non-evolving galaxy models for B > 2A (eg. 
Tyson 1988). Simulations based on star-formation rates that assume yields in metallicity from 
local observations, predict that the amount of metals (and hence dust, assuming a hxed fraction 
of metals condense into grains at a constant rate) produced from such a population will be smaller 
than local estimates by an order of magnitude (eg. White & Frenk 1991). If the only source 
of IGM dust was from these ‘low-mass’ galaxies, then the total optical depth to z>2 would be 
insignihcant, and effects on the background universe would be minimal. 


6 Conclusions 

In this paper, we have shown that dust is more effective at obscuring background sources when 
diffuse or extended. We find that obscuration of background sources by a given dust distribution 
with optical depth tb will be most effective when tb < '^/13, where (3 is the cumulative luminosity 
function slope of the sources. 

We have explored the effects of diffuse dust from, first, galaxy clusters and second, from 
a hypothesised uniform IGM component. By assuming different radial dust density prohles in 
a typical rich cluster, we have predicted the optical depth and background source dehcit as a 
function of projected cluster radius. These predictions can be compared with future observations 
to constrain the properties of intracluster dust. Our predicted optical depth measures (ts< 0.3) 
satisfy the above criterion {tb < 1//5) for background luminous QSOs and galaxies. Existing 
studies claiming anticorrelations in the distribution of QSOs with foreground clusters down to 
scales ~ F are consistent with a dust density prohle that follows the galaxy distribution. 

As a further illustration, we have explored the effects of a diffuse IGM dust component with 
cosmic mass density equal to that observed in local galaxies. Assuming this density is constant 
on comoving scales to z = 2, we hnd a dehcit in background sources about three times greater 
than that predicted assuming dust in normal galaxies alone. 

The ‘dihuseness’ of the dust is the key parameter which we claim determines the ehectiveness 
of obscuration of the background universe. Although such dust distributions may be difficult to 
detect, we must not neglect their possible presence. Further studies of spatial dust distributions, 
preferably via the counts and colours of background sources will be essential in conhrming our 
predictions. 
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